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1. Introduction 



Spectral functions of mesonic operators play an important role in finite temperature 
QCD. Many experimental results in high energy heavy ion collisions ( e.g. low invariant 
mass dilepton enhancement, anomalous J/ijj suppression etc.) can be understood in 
terms of medium modifications of meson spectral functions P]. 

For a long time it has been believed that lattice QCD is suitable only for calculation 
of static quantities at finite temperatures, such as the transition temperature, equation of 
state, screening lengths etc. However, it was shown by Asakawa, Hatsuda and Nakahara 
that using the Maximum Entropy Method one can in principle reconstruct also meson 
spectral functions. The method was successfully apphed at zero temperature [3 13 El- 
and also at finite temperature [3 El 13 IHl El 113 HI] • Though systematic uncertainties 
t speaker at the conference 
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in the spectral function calculated on lattice are not yet completely understood, it was 
shown in Ref. j3] that precise determination of the imaginary time correlator can alone 
provide constraints on the spectral function at finite temperature. 

In this contribution we are going to review recent results on meson correlators and 
spectral functions in finite temperature QCD. 



2. Meson spectral functions 

The spectral function for an operator is directly related to the Fourier transforms of real 
time correlators of the operator. Defining the real time propagators 

D<{t) = (6(0)O(t)), D>{t) = (O(t)O(O)), (1) 

for a suitable operator O, the spectral function a{uj) is defined as 

D>iu)-D<{u) 

^(^) = ^ ' (2) 

where D^''^{uj) are the Fourier transforms of the correlators in Eq. (Q). The spectral 
function encapsulates the properties of the channel corresponding to operator O: a 
stable mesonic state will show up as a 5 function and an unstable state will be present 
as a broad (Breit-Wigner) peak. Also, other quantities of experimental interest can 
be directly extracted from spectral functions: for example, the rate of emission of 
thermal dileptons can be expressed in terms of the spectral function of the vector current 

dW 5aL 1 



em 



dpod^p - 27n^ pliePo/T _ ^fviPo,P)- (3) 

In lattice studies of finite temperature systems in equilibrium, on the other hand, 
one calculates the imaginary time or Matsubara correlator 

G{r) = (O(r)6(0)), (4) 

where the Euclidean time r G (0,/3 = 1/T). The imaginary time correlators is an 
analytic continuation of the real time correlator G{t) = D^{—iT). This and the 
periodicity in r lead to an integral representation for G{t) in terms of the spectral 
function, 

POO 

G{t,T) = / duja{uj,T)K{T,uj,T) (5) 

).2^!iMi^. (6) 

smh(co'/ Zi ) 

Here we consider correlators of point meson operators qTq, with F = 1, 75, '~ffj,,'~f^'y5 
for scalar, pseudo-scalar, vector and axial-vector channels. We will also restrict the 
discussion to the case of zero spatial momentum. 

Equation Q is valid in the continuum. It is not obvious that it holds also for 
the correlators calculated on the lattice. However, it was shown in that an integral 
representation like in Eq. (jHl) can be also defined for the lattice correlator in the limit of 
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the free field theory, with the only modification being that the upper integration limit 
becomes finite Wmax — 4a~^ ^Hl- The effect of the discretization is contained solely in 
the spectral function 

In extracting a{uj) from the lattice measurements of G{t) using Eq. the obvious 
problem is that we have to reconstruct several hundred degrees of freedom (needed for a 
reasonable discretization of the integral in Eq. (jSj) with (9(10) data point on G{t). This 
ill-posed problem can be solved using the Maximum Entropy Method (MEM) [2]. In 
the Maximum Entropy Method on searches for spectral function which maximizes the 
conditional probability exp(— x^/2 + aS'), where is the standard likelihood function 
for the data and 



is the Shannon- Jaynes entropy which parametrizes prior knowledges such as positivity 
of the spectral functions and informations about the asymptotic form of the spectral 
function. Finally a is a real and positive parameter, which is integrated out to obtain 
the most probable spectral function. The Shannon- Jaynes entropy depends on one real 
function m{u) called the default model. 

3. Numerical results for the correlators and spectral functions 

Here we will mostly discuss finite temperature spectral functions for heavy quarkonia. 
The behavior of heavy quarkonia at finite temperature has been of great interest ever 
since the suggestion of Matsui and Satz that the disappearance of J/ip can act as a 
signal of onset of deconfinement, since the screening in a quark-gluon plasma will lead 
to the dissolution of charmonium bound states IMJ. The problem of the J/ip suppression 
was addressed using potential models [HI [13 [IHI, and resulted in estimates that the 
J/ip will dissolve around l.lTc [inilinj. More recent calculations, using lattice results for 
internal energy of a static qq pair [TJ] as the potential, lead to a much higher dissociation 
temperature > l.ST^ for the J/ip [THIIT^ . 

Direct lattice investigations of heavy quarkonia at finite temperatures, by studying 
the Matsubara correlators and extracting the spectral funtion from them using MEM, 
were performed only recently. Calculations have been done using anisotropic [HI Elj 
as well as fine isotropic [TT] lattices in quenched approximation (i.e., for a purely 
gluonic plasma). A check of the cutoff dependence of the results by using several lattice 
spacings, which turns out to be quite important, has been done so far only in Ref. [TT]. 
The extracted spectral functions at zero temperature reproduced the properties of the 
ground state mesons, but revealed structures (lattice artifacts) at higher energies whose 
positions scale like [31[11[TT]. Due to these structures, the high energy part of the 
spectral functions in the interacting theory look very different from that in the free 
theory. 

We first discuss the spectral functions at low temperatures where the large extent 
of the imaginary time direction allows a reliable extraction of the spectral function from 
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Figure 1. The charmonia spectral functions below deconfinement for two different 
lattice spacings a"! = 4MGeV (left) and a''^ = 9.72GeV (right) [TT] . 
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Figure 2. Ratios G/Grecon for the pseudo-scalar and the vector correlators (left) and 
for the scalar and the axial- vector correlators (right). 

the Matsubara correlators. We have used the continuum free meson spectral function 
of the form mou'^ as a default model. It was checked that the result does not depend on 
the choice of the default model ^T]. In Fig. ^we show charmonium spectral functions 
below deconfinement at two different lattice spacings, corresponding to temperatures of 
0.75Tc and 0.9Tc [TT]. The ground state peaks in the pseudo-scalar and vector channels 
correspond to r]c and J/ip states, while those in the scalar and axial- vector channels 
correspond to Xco and Xci states, respectively. While the ground state peaks reproduce 
the expected masses the positions of the other peaks scales roughly as 1/a, i.e. these 
peaks are lattice artifacts. 

Now let us discuss what happens at finite temperature. In this case the 
reconstruction of the spectral function becomes more difficult as not only the number 
of available time-slices is reduced but also the extent of the imaginary time direction 
(= 1/T) may be quite small (see discussion in Ref. P on this point). While the 
former limitation can be relaxed using larger anisotropies, the latter will always remain. 
Therefore it is useful to inspect the temperature dependence of the correlators first. 
The temperature dependence of the correlators is due both to the explicit temperature 
dependence of the corresponding spectral functions and to the temperature dependence 
of the integration kernel in Eq. (jH)). To factor out the trivial temperature dependence 
due to the kerneLwe construct the model correlators 



POO 

Grec{T,T) = / duja{uj,T* 
Jo 



)K{t,u;,T), 



(8) 



with T* being some temperature below deconfinement. If the spectral function has 
no significant temperature dependence also above deconfinement, we would expect 
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Figure 3. Charmonia spectral functions for the vector (left) and scalar (right) channels 
above Tc. mj. 

G/Grec — 1 also at temperatures above Tc. In Fig. |2l we show G/Grec for different 
channels at 1.5Tc. As one can see this ratio is not very far from unity for pseudo-scalar 
and vector channels, while for the scalar and axial-vector channel large changes are 
visible. In fact, for the scalar and axial- vector channels large changes are seen already 
at l.lTc [TTj. This suggests that the IS state charmonia (r/c and J/ip) may survive 
deconfinement with little modifications of their properties, while IP state charmonia (xco 
and Xci) are dissolved or strongly modified by the medium right after deconfinement. 
Notice that Fig. |21 shows the ratio G/Grec for two different lattice spacings, with very 
little lattice spacing dependence observed in the results. 

To study the problem in more detail one needs to extract the spectral functions 
from the measured correlators. Because of the limited extent of the time direction 
this becomes quite difficuilt at high temperatures, as the result gets sensitive to the 
default model. We have used the high energy part of spectral function below Tc as 
part of the default model. This is justified as the high energy part of the extracted 
spectral function is dominated by lattice artifacts, and also, the effects of moderate 
temperatures should not be significant at very high energies. In Fig. |S1 we show the 
spectral functions extracted for the vector and scalar channels, at temperatures of 1.5 
- 3 Tc. The figure confirms the expectation based on the analysis of the correlator; 
the J/ip is seen to exist up to temperatures 2.25Tc , while the Xco state is dissolved 
already at LIT^. The situation is similar for the pseudo-scalar and axial vector spectral 
functions. The results for the survival of the IS" states till quite high temperatures 
are confirmed by calculations done on anisotropic lattices [HI CD] , though there is some 
controversy about the temperature where J/ip dissolves and the way it disappears. In 
Ref. it was found that the 15* charmonia abruptly disappear at 1.7Tc, while the 
spectral functions shown above suggest a gradual dissolution with resonance structure 
surviving till 2.25Tc. 



4. Conclusions 



We have discussed lattice calculations of the charmonium spectral functions at finite 
temperature. Despite the uncertainties present in the current analysis we can conclude 
that J/ip and rjc survive in the gluonic plasma up to quite high temperatures, possibly to 
temperatures as high as 2.25Tc. Moreover, the properties of the IS charmonia seems to 
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be only very mildly affected by the deconfinement transition. The IP states charmonia 
on the other hand are dissolved right after deconfinement. This conclusion is supported 
both by the explicit calculations of the spectral functions using MEM and by the analysis 
of the correlators which is much more robust. The spectral functions calculated on the 
lattice show strong lattice artifacts at high energies. However the properties of the 
first peak, which corresponds to the physical state, are not strongly affected by these 
artifacts. This conclusion is also supported by direct analysis of meson correlators at 
two different lattice spacings. 

Finally we note that meson spectral functions have also been analyzed for the case of 
light quarks El IE] • Contrary to the case for the heavy quarks, the spectral functions 
for light mesonic observables show quite strong temperature dependence. However, 
also here one sees resonance like structures at low energies and no sign of free quark 
propagation. The lattice artifacts at high energy are present also in this case but were 
not studied in detail yet. This would be absolutely necessary for drawing physical 
conclusions with confidence. 

Acknowledgments 

This work has been partly supported by the U.S. Department of Energy under contract 
number DE-AC02-98CH10886. P.P is Goldhaber and RIKEN Fellow. 

[1] Rapp R, Wambach J 2000, Adv. Nucl. Phys. 25 1; Rapp R, Grandchamp L, hep-ph/0305143 

[2] Y. Nakahara, M. Asakawa and T. Hatsuda, Phys. Rev. D 60 (1999) 091503; Prog. Part. Nucl. 
Phys. 46 (2001) 459 

[3] T. Yamazaki et al. [CP-PACS CoUaboration] , Phys. Rev. D 65 (2002) 014501 

[4] T. Blum and P. Petreczky, arXiv:hep-lat/0408045 

[5] F. Karsch, E. Laermann, P. Petreczky, S. Stickan and I. Wetzorke, Phys. Lett. B 530 (2002) 147 

[6] F. Karsch, S. Datta, E. Laermann, P. Petreczky, S. Stickan and L Wetzorke, Nucl. Phys. A 715 
(2003) 701 

[7] M. Asakawa, T. Hatsuda and Y. Nakahara, Nucl. Phys. A 715 (2003) 863 

[8] P. Petreczky, J. Phys. G 30 (2004) S431 

[9] T. Umeda, K. Nomura and H. Matsufuru, arXiv:hep-lat/0211003 

[10] M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 92 (2004) 012001 

[11] S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, Phys. Rev. D 69 (2004) 094507 

[12] E. Braaten, R. D. Pisarski and T. C. Yuan, Phys. Rev. Lett. 64 (1990) 2242 

[13] F. Karsch, E. Laermann, P. Petreczky and S. Stickan, Phys. Rev. D 68 (2003) 014504 

[14] T. Matsui and H. Satz, Phys. Lett. B 178 (1986) 416. 

[15] F. Karsch, M.T. Mehr, H. Satz, Z. Phys. C 37 (1988) 617 

[16] S. Digal, P. Petreczky and H. Satz, Phys. Rev. D 64 (2001) 094015, Phys. Lett. B 514 (2001) 57 

[17] O. Kaczmarek, F. Karsch, P. Petreczky, F. Zantow, Nucl. Phys. B (Proc. Suppl.) 129 (2004) 560 

[18] E. V. Shuryak and I. Zahed, arXiv:hep-ph/0403127. 

[19] C. Y. Wong, arXiv:hep-ph/0408020 



